Angular momentum content of the p-meson in lattice QCD 
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The variational method allows one to study the mixing of interpolators with different chiral 
transformation properties in the non-perturbatively determined physical state. It is then possible to 
define and calculate in a gauge- invariant manner the chiral as well as the partial wave content of the 
quark-antiquark component of a meson in the infrared, where mass is generated. Using a unitary 
transformation from the chiral basis to the LSJ basis one may extract a partial wave content of a 
meson. We present results for the ground state of the p-meson using quenched simulations as well 
as simulations with n/ = 2 dynamical quarks, all for lattice spacings close to 0.15 fm. We point out 
that these results indicate a simple 3 5i-wave composition of the p-meson in the infrared, like in the 
SU(6) flavor-spin quark model. 

PACS numbers: 11.15.Ha, 12.38.Gc 



1. Motivation. The composition of hadronic states in 
quantum field theory is a subtle issue. Whereas in non- 
relativistic approaches the notion of a wave function and 
a complete basis of states is well-defined, in quantum field 
theory beyond the ground state there is no well-defined 
single hadron, and the state is always a scattering state 
with superposition of many particle components. A given 
hadron interpolator couples in principle to all states with 
its quantum numbers. 

Lattice studies of hadrons are typically limited to spec- 
troscopy, static observables like magnetic moments, ax- 
ial couplings, etc., as well as dynamical observables such 
as form factors, parton distributions, etc. A challeng- 
ing task is to reveal the composition of hadrons, i.e., to 
understand the hadron structure in ab-initio QCD cal- 
culations. Of course, in principle any hadron contains 
a large amount of different Fock components and it is 
hardly possible to reconstruct on the lattice the complete 
hadron wave function. Nevertheless, the analyses of the 
phenomenological data and modeling of hadrons suggest 
that typically only a few components are the dominant 
ones. To understand in ab-initio calculations the struc- 
ture of the leading component in the infrared, i.e., at 
the scale where hadron mass is generated, would be an 
important step to improve the insight into hadron struc- 
ture. 

There is a tool to study the hadron wave function 
on the lattice - the variational method [l], B (for a 
recent, more complete set of references see [3j). The 
correlation function of a hadron interpolator, which is 
built from a linear combination of many interpolating 
operators with correct quantum numbers will single out 
an optimal signal associated with an exponential de- 
cay for large Euclidean time distance. This combina- 
tion defines the "physical" ground state hadron. In 
the variational method one chooses a set of interpola- 
tors {Oi, O2, On} that potentially couple to a given 
hadron and computes the correlation matrix 
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If the set is complete enough the eigenvalues of a general- 
ized eigenvalue problem may be related to the eigenstates 
of the Hamiltonian. Then, if one is interested in recon- 
structing the leading Fock component of a hadron in the 
infrared, e.g., the qq component of a meson, one needs 
interpolators that allow one to define uniquely such a 
component. 

Let us consider the chiral limit of two- flavor mesons. 
All possible quark-antiquark interpolators for non-exotic 
mesons can be classified according to representations of 
the SU(2) L xSU(2) R and U(1) A groups This basis 

is complete for a quark-antiquark system provided that 
there is no explicit excitation of the gluonic field with 
non- vacuum quantum numbers. Using interpolators that 
belong to this basis allows one to study the chiral sym- 
metry breaking aspects in a meson. Diagonalizing the 
cross-correlation matrix one can reconstruct a decompo- 
sition of a given meson in terms of different representa- 
tions of the chiral group. Chiral symmetry breaking in 
the infrared would mean that the meson wave function 
has components with different transformation properties 
with respect to SU(2) L x SU{2) R and U(1) A - 

It is also of interest to reconstruct a composition of a 
meson in terms of the 2S+1 Lj basis, where J = L + S 
are the standard angular momenta. Such a decompo- 
sition provides a bridge to the language of the quark 
model. A priori it is clear that in the heavy quark 
limit the non-relativistic language of the quark model is 
nearly adequate. There are well known achievements of 
the quark model, such as the SU (6) flavor-spin symme- 
try, which is related to rather large masses of the con- 
stituent quarks. These large masses are far from the 
tiny masses of the current quarks and it is unclear what 
happens for light quarks. While there is a kind of un- 
derstanding that successes of the quark model are re- 
lated to chiral symmetry breaking and that these large 
masses of constituent quarks arise from the coupling of 
current quarks to the quark condensate at low momenta, 
a model-independent and gauge-invariant view on this 
problem is absent. Then, it would be intriguing to see a 
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decomposition of the leading quark-antiquark Fock com- 
ponent of a meson in terms of the 2S+1 Lj basis in the 
infrared. 

There is a possibility to establish the angular mo- 
mentum decomposition of the leading quark-antiquark 
component. Both the chiral basis and the 2S+1 Lj ba- 
sis are complete for a two-particle system. There ex- 
ists a unitary transformation from the chiral basis to the 
standard 2S+1 Lj basis in the center-of-momentum sys- 
tem [f|. Each of the states in the relativistic chiral ba- 
sis can be uniquely represented in terms of the allowed 
2S+1 Lj states. Diagonalizing the cross-correlation ma- 
trix in terms of the interpolators with carefully chosen 
chiral transformation properties, and using this unitary 
transformation to the 2S+1 Lj basis one can reconstruct 
a partial wave decomposition of the leading Fock compo- 
nent of a meson. 

The method is general and can be applied to any me- 
son. Here we use as an example the p-meson. If chi- 
ral symmetry is unbroken, then there are two possible 
states in the chiral basis with the p-meson quantum num- 
bers, |(0, 1)+(1, 0); 1 1— ) and |(l/2, l/2) 6 ; 1 1~ ), where 
(0, 1) + (1,0) and (1/2, l/2)& specify two different repre- 
sentations of SU(2)l x SU(2)r that are compatible with 
the p-meson quantum numbers /, J PC = 1,1 (a sub- 
script b specifies one of the two different representations 
(1/2, l/2) and (1/2, l/2) 6 0,(1). These two representa- 
tions are a complete set for particles with the quantum 
numbers of the p-meson. A unitary transformation from 
the chiral basis to the 2S+1 Lj basis takes the form 

1(0,1) + (1,0); 11") = y||l; 3 ^) + y||l; 3 A), 
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Hence, if we choose two different interpolators with 
the same spatial content (i.e., with same smearing width 
of the quarks; its size defines the scale where we probe 
the hadron) and transformation properties according to 
|(0,1) + (1,0);1 1— ) and |(l/2, l/2) b ; 1 l"), we will be 
able to reconstruct the angular momentum content of the 
qq component. Such interpolators are well-known - the 
quark bilinears qj' l rq and qa 0l rq 7]. It has been estab- 
lished in lattice simulations that the p-meson couples to 
both (see, e.g., [1,0, E3] for quenched and 11 1 for dynam- 
ical simulations). A key property is that these two in- 
terpolators have radically different chiral transformation 
properties. The former one transforms as (0, 1) + (1, 0) 
while the latter one belongs to the (1/2, 1/2);, representa- 
tion of SU (2) l x SU (2)_r. Consequently, by diagonalizing 
the cross-correlation matrix with the effrq and qa 0l Tq 
interpolators we can reconstruct a decomposition of a 
given eigenstate (i.e., of the ground state p-meson and 
its excitations) in terms of the possible partial waves in 



the infrared. 

2. Analysis of lattice correlators. For completeness 
we briefly summarize the features of the variational anal- 
ysis that allow the decomposition of the ground state 
properties in the p-channel. 

The normalized physical states \n) propagate in time 
with 



(n(t)\m(0)) =6„ 



(3) 



The interpolating (lattice) operators Oi{t) are projected 
to vanishing spatial momentum and are usually not nor- 
malized. We compute the correlation function 
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with the coefficients giving the overlap of the lattice op- 
erator with the physical state, 
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For interpolating operators Oj spanning an orthogonal 
basis these values would indeed constitute the wave func- 
tion of state |n(0)) in that basis. 

We assume that the correlation matrix §Z§ can be ap- 
proximated by a finite sum over N states and denote this 
approximation by CCt)ii 



It can be shown 0, | (for a recent discussion see 0, 
12j ) that the generalized eigenvalue problem 

C(t) t3 uf =\W{t,t )d(tohv$ ,) (6) 

allows to recover the correct eigenvalues and eigenvector 
within some approximation. One finds 

A(")(i,t ) = e- £( " > ( t - t «)(l + o(e- A£( " , ( t ^))) , (7) 

where AE^ may be as small as the distance to the next 
nearby energy level. In 0] it was pointed out that in an 
interval to < t < 2 to these contributions are suppressed 
and leading terms even have AE^ equal to the distance 
to the first neglected energy level En+i- At to all eigen- 
values are 1 and the eigenvectors are arbitrary. 

The eigenvectors come out orthogonal (dual) to the 
original wave functions a'"', 
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and approximate the correct ones. Here o- m > is a nor- 
malization which we get rid off below. 
We define a sum of lattice operators 



? (*)=£ u (n)* . 



(9) 
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and find 



N 
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Therefore 



jjWtlQ) = c (")*|„) 
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The eigenvector coefficients are related to the composi- 
tion of the eigenstate in terms of the interpolating oper- 
ators. The original values af 1 are then 

of = (OlOiln) = -1^.(010* rjWt| Q) 

AT 



Jn) 



They would agree with uf^ if the interpolators were or- 
thogonal, which they are usually not. 

However, with and (JSJ) we find for the large t be- 
havior (summation convention) 
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and for the ratio for large t 
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Assuming asymptotically leading exponential behavior 
this allows to read off \a\ \ in the asymptotic region. 
The ratio 



(15) 



tells us how much the interpolating operator Oi con- 
tributes to the eigenstates \n) and |m). This can be 
used to discuss the ratio of decay constants of various 
excitations as done in 0, [T^, [l3| • 

If we are interested in ratios of couplings of the different 
lattice operators to the physical states, we may utilize 



C(t)ijU 
C{t) k3 u 
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This ratio tells us how much different interpolating oper- 
ators contribute to the eigenstate \n). This can be used 
to discuss contributions of, e.g., different representations 
of the vector meson channel. 



TABLE I: Specification of the data used here; for the gauge 
coupling only the leading value /3lw is given, mo denotes the 
bare mass parameter of the Cl-action. Further details like the 
determination of the lattice spacing and the 7T- and p-masses 
are found in [l5l . Il6l ] . For the quenched case and ensemble A 
we used 100 configurations, for sets B and C we analyzed 200 
configurations each. The lattice size is 16 3 x 32. 



Data 



Pi 



a mo 



[fn 



r [MeV] m p [MeV] 



Quenched 7.90 0.04-0.20 0.1480(10) 475-1053 912-1251 



dyn. 
dyn. 
dyn. 



A 4.70 
B 4.65 
C 4.58 



-0.050 
-0.060 
-0.077 



0.1507(17) 
0.1500(12) 
0.1440(12) 



526(7) 922(17) 
469(4) 897(13) 
318(5) 810(28) 



3. Lattice simulation and results. In a series of pa- 
pers we have studied the hadron spectrum derived in 
the quenched case [IF 
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|15| | as well as for dynamical 
fermions [16j. The gauge field action was the Liischer- 
Weisz action and the fermions were simulated with the 



so-called chirally improved (CI) Dirac operator [171. Il8j. 

In these analyses the variational method was used; 
the hadron interpolators were built from smeared quark 
sources. Here we refer only to a subset of these results, 
namely the p-channel ( J PC — 1 ) and restrict ourselves 
to the four interpolators 

Ol = Unl l d n , 02 = Uwl'dw , (17) 
03 = UnJ^dn , 4 = U^fdv, . (18) 

Here Y is one of the spatial Dirac matrices, 74 is the 7- 
matrix in (Euclidean) time direction, and the subscripts 
n and w (for narrow and wide) denote the two smearing 
widths, 0.25 fm and 0.41 fm, respectively [151 ]. 

All results have been obtained on lattices of size 16 3 x 
32 with lattice spacing a close to 0.15 fm (see Table HJ. 
The runs A, B and C are for n/ = 2 mass degenerate 
Cl-fermions. 

As detailed in [lil, [l6[ it is possible to identify the 
ground state and the first excited state. In that study 
we also presented the corresponding eigenvector compo- 
nents (u^ in our present notation) which were stable 
over several time slices. Continuing that analyses we 
now utilize Eq. (|16[) and determine the ratios ai/a k for 
the ground state (we have dropped the superscript refer- 
ence to the physical state), in particular 01/03 and 02/04 
which give information on the relative contributions of 
operators with 7* compared to 7*7*. 

We chose the relatively small value to = 1. Larger 
values lead to larger errors; however, we found very nice 
plateau behavior over a range of propagation distance 
in Euclidean time t. In Fig. [T] we show the surprising 
stability of the plateau of the ratio 01/03 for the three 
data sets obtained for dynamical fermions. We average 
the values in the range 3 < At < 7 and show the resulting 
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FIG. 1: As an example we show the ratio of the contributions 
of Oi vs. O3 (the narrow-narrow smeared operators) for en- 
sembles A, B and C. All four operators have been included 
in the correlation matrix. We do not show the points for 
At < to = 1 and above At — 10, where approaching the time- 
symmetry points the results become statistically unstable. 
The error bars have been determined with single-elimination 
jack-knife. 
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FIG. 2: We show the plateau mean values of the ratios 0,1/0,3 
and 02/ '0,4 for the quenched (circles) and dynamical (filled 
squares) data sets considered in Table [I] For the dynami- 
cal data we also have partially quenched results, i.e., with 
m-vaience > rrisea (not shown here), which approach the purely 
quenched values for pion masses above «0.7 GeV. The error 
bars have been determined with single-elimination jack-knife. 



ratio in Fig. [2] for both, the quenched and the dynamical 
data. We find different behavior towards the chiral limit, 
indicating the effect of the fermion sea. 

Extrapolation towards the physical point indicates a 
value close to 1.35, i.e., compatible with \/2 within the 
errors. Consequently chiral symmetry is broken such 
that the qq component of the p-meson is a superposi- 
tion of the (0, 1) + (1,0) and the (^, i) h representations 
with the ratio close to Inverting the unitary ma- 

trix in ([2|) we then conclude that the physical p-meson 
in the infrared is almost purely a Si state. This gives 
a model-independent and gauge invariant explanation of 



the success of the SU(6) flavor-spin symmetry for the p- 
meson. 

This result holds in the infrared where the mass is gen- 
erated. We probe the chiral and partial wave content 
of the p wave function at a scale fixed by the smearing 
size. At higher resolution (Q 2 — ► 00) the /3-function ap- 
proaches zero, the pseudotensor current decouples from 
the p-meson and the partial wave decomposition should 
be determined by the (0, 1) + (1, 0) representation. Such 
a tendency is indicated in Fig. [21 It is also consistent 
with the ratio ~ 1.6 obtained in dynamical calculations 
with point interpolators (i.e., a scale fixed by a — 0.114 
fin) in [ll|. 
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